Viscoelastic Depinning of Driven Systems: Mean-Field Plastic Scallops 
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We have investigated the mean field dynamics of an overdamped viscoelastic medium driven 
through quenched disorder. The model allows for the coexistence of pinned and sliding regions and 
can exhibit continuous elastic depinning or first order hysteretic depinning. Numerical simulations 
indicate mean field instabilities that correspond to macroscopic stick-slip events and lead to pre- 
mature switching. The model describes the elastic and plastic dynamics of driven vortex arrays in 
superconductors and other extended disordered systems. 



Extended condensed matter systems driven over 
quenched disorder exhibit a very complex dynamics, in- 
cluding nonequilibrium phase transitions and history de- 
pendence. Such systems include vortex arrays in type-11 
superconductors charge density waves in anisotropic 
conductors H,^, and many others. Closely related be- 
havior also arises in friction and lubrication Q, where a 
surface or monolayer is brought in contact with another 
solid surface and forced to slide relative to it. 

Most of the theoretical work to date has focused on 
the dissipative dynamics of driven elastic media that are 
distorted by disorder, but cannot tear. At zero temper- 
ature such systems exhibit a sharp depinning transition 
from a pinned to a sliding state The transition, 

first studied in the context of charge density waves, is 
continuous, with universal critical behavior. The sliding 
state is unique and there is no hysteresis or history de- 
pendence More recent work, while still focusing on 
elastic media, has shown that the dynamics is quite rich 
well into the uniformly sliding state [f7|-pd| . 

On the other hand, experiments |jl^,|l3| and simula- 
tions show that the elastic medium model is in- 
adequate for many physical systems with strong disorder 
that upon depinning exhibit a spatially inhomogeneous 
plastic response, without long wavelength elastic restor- 
ing forces. In this plastic flow regime, topological de- 
fects proliferate and the system is broken up in fluid-like 
regions flowing around pinned solid regions. Not much 
progress has been made in describing this behavior ana- 
lytically. The wealth of experimental work on driven vor- 
tex arrays clearly indicates that, in most of the field and 
temperature region of interest, the current-driven vortex 
dynamics is strongly history dependent, with long-term 
memory and switching as the system explores a variety 
of nonequilibrium sliding states ||l^,|l3| . 

In this paper we describe a coarse-grained model for 
the dynamics of a driven medium that allows for spatially 
inhomogeneous response, with the coexistence of moving 
and pinned regions. The model is inspired by the well- 
known phenomenology of viscoelasticity in dense fluids 
p6[. The elastic couplings between the local displace- 



ments are replaced by couplings that are nonlocal in time 
and allow for clastic restoring forces to turn into dissipa- 
tive fluid flow on short time scales. The model yields 
elastic depinning in one limit; as the parameters are 
varied, it incorporates continuous depinning, hysteretic 
plastic depinning and eventually viscous flow, allowing 
the crossovers between these regimes (such as those, ob- 
served in vortex arrays ||l2|,|l5|) to be studied in detail. 
For a wide range of parameter values the depinning tran- 
sition is first order, with velocity hysteresis (switching.) 
The nonlinear velocity-force characteristic can be evalu- 
ated analytically in mean flcld for various types of pin- 
ning forces, under the assumption of constant mean field 
velocity. Numerical simulations confirm the inhomoge- 
neous nature of the dynamics, with pinning and tearing 
(coexisting moving and pinned degrees of freedom.) In 
addition, the mean velocity near depinning fluctuates, 
due to macroscopic stick-slip type events. These events 
appear to only mildly violate the uniform mean-velocity 
assumption but directly lead to switching from one veloc- 
ity branch to another before the first branch terminates 
(premature switching.) Models that account for switch- 
ing in charge density waves and are in spirit similar to 
ours have been proposed and studied by Strogatz and 
collaborators [0. In such models, plasticity is modeled 
by a non-convex elastic potential, in contrast with the 
velocity convolutions studied here. A model similar to 
ours has also been proposed for crack propagation ||l8| . 

The model: a driven viscoelastic medium. To motivate 
our model, we first recall the generic model of driven elas- 
tic media |^] discussed extensively in the literature, where 
the long-wavelength dynamics is described in terms of a 
coarse-grained displacement field, u{r,t). The displace- 
ment fields represent deformations of regions pinned col- 
lectively by disorder (e.g., a Larkin domains) and are 
coupled by convex elastic interactions. No topological 
defects are allowed. Considering, for simplicity, the over- 
damped dynamics of a scalar field (the model is easily 
extended to the more general case) and modeling the 
displacement field on lattice sites, u{r,t) — > Ui{t), the 
equation of motion for the local field Ui (measured in the 
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laboratory frame p9|) at site i is 
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fiij {uj -Ui) + F + Fi{ui), 



(1) 



where the summation is restricted to nearest neighbor 
pairs and 70 is the friction. If all the nearest-neighbor 
elastic couplings, fiij > 0, are equal, the first term on 
the right hand side of Eq. (Q) is the discrete Lapla- 
cian in d dimensions. The second term is the homoge- 
neous driving force, F, and Fi(ui) denotes the pinning 
force arising from a quenched random potential, Vi{ui)^ 
Fi{ui) = -dVi/dui = hif{ui - (3i), with f{u) a periodic 
function with period 1 and Pi random phases uniformly 
distributed in [0, 1]. The hi are chosen independently at 
each site from a distribution p{h). One of the quantities 
of interest is the average velocity of the driven medium, 
v{F) — N~'^'Y^Ui. For an elastic medium there is a 
unique stationary sliding state for F > Fc, with critical 
behavior v(F) ^ {F — Fc)'^ and no hysteresis at the 
transition 

We now modify the elastic interactions in Eq. (|^) to 
allow for local tearing of the medium. Inspired by stan- 
dard models of viscoelasticity, we replace the elastic in- 
teraction by couplings to the local velocity field, Vi = iii, 
that are nonlocal in time. Our model equation for the 
overdamped dynamics of a driven viscoelastic medium is 
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where the viscous couplings Cij{s) have finite first mo- 
ments, /p dsCij{s) = rjij < 00 and Cij(O) = /iy. Such 
nonlocal couplings to velocity are of course not present 
at the microscopic level, but are generated generically 
upon coarse-graining ||l6| , ^ . Eq. (g) is a coarse-grained 
model for the dynamics of a driven disordered medium 
that allows for slip or friction of the interacting Larkin 
domains relative to each other. 

A simple, yet successful, model of viscoelasticity due 
to Maxwell is obtained when the memory kernels are as- 
sumed to be uniform in space and to decay exponentially 
in time, according to Cy (t) = /ze" 
Maxwell relaxation time. For r 

reduces to Eq. for a driven elastic medium. For 
T — > and 77 fixed, the first term on the right hand side 
of Eq. (||) can be approximated as ilj^^ij) b^j{t) ~ 
which represents viscous forces coupling the local fluid 
velocity at different spatial points. In this limit, Eq. (j^) 
describes an overdamped lattice-fluid of viscosity rj. We 
propose Eq. (^) as a simple, yet realistic model for a 
driven disordered system that exhibits spatially inhomo- 
geneous plastic response. 
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Mean Field Approximation. As for the driven elastic 
media, substantial analytical progress in two or three di- 
mensions is presumably only possible via perturbation 
theory or by a functional renormalization group treat- 
ment p^ . An alternative approach that has provided 
valuable insight for a driven elastic medium is mean field 
theory (MET), first discussed by D. S. Fisher |]. MET 
is formally valid in the limit of infinite-range interaction, 
with J2j C-ij — NC{t) held fixed. The equation of motion 
for the displacement at each site is then given by 

7oUi= / ydsC{t^ s)[v{s) ^u,{s)] + F + F,{ui), (3) 



where the mean field is given by u{t) — iV^^ Si^i Ui(Oj 
and v(t) = u{t). 

If the memory kernel C{t) is chosen to be of the 
Maxwell form, it is thcnpossible to transform the integro- 
differential equation to a second-order differential 
equation, given by 



rit -I- 7(77, T, h; Ui)u = F + Fi{ui) +r]v, 



(4) 



with "f{ri,T,Ui;h) = 1 + 7] — '^^^ effective friction. 
We have scaled Eq. (Q) by letting r tHq, t thQ, 
V v/lo, F F/iloho) and h h/{'yoho), where ho 
is the characteristic scale of the distribution p{h). With 
this change of variables, the model is now characterized 
by two parameters, 77 and r, and the shape of p{h). The 
MF equation for our viscoelastic model closely resembles 
the MF equation for a driven massive elastic medium, 
with T playing the role of the mass. The most important 
difference is that in the massive elastic medium the MF 
term rjv is replaced by fiu. As a result, the MET of a 
driven massive elastic medium even with constant v con- 
tains three degrees of freedom (as opposed to the two of 
our problem) and the dynamics of a single Ui is chaotic 

0,0 ■ 

We are first interested here in steadily sliding solutions 
of the MF model, Eq. (|). It is natural to look for peri- 
odic solutions Up{t; h) of period T{h), {Jq'^'^^ dtup{t; h) = 
1) that may set in after an initial transient {t ^ T(h)). 
Such solutions need not be unique. Guided by a large 
body of previous work on driven elastic media, we focus 
on the MET for the case of a pinning potential with cusp- 
like singularities, which better captures the physics of 
the corresponding finite-dimensional model ||23| . An ex- 
plicit solution of Eq. (^) can be obtained for the scalloped 
parabolic potential, V{u) = {h/2){u'^ — m + 1/4). In this 
case Eq. is linear and its general solution is Up{t; h) = 
Ci exp(-Aii) + C2 exp(-A2t) + 1/2 + {rjv + F)/h, with 
Ai,2 = (l + T] + Th± + ?/ + t/7)2 - 4t/i) /(2r). For 
each fixed value of h, we obtain an implicit equation for 
the period, T{h), jfu + F = G(T] rj, r, ft,), with 
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The solution of Eq. , together with the self-consistency 
constraint v =< [T{h)]^^ >h, determines the drift veloc- 
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ity as a function of driving force, F. When T{hi) 
the Ui is pinned. 

Figure 1 shows the analytical solution for the mean ve- 
locity as a function of driving force for p{h) — exp(— ft,). 
The depinning occurs at F — for all distributions of 
pinning strengths, p{h), with support not bounded from 
below by a positive /imin- For small rj and r, correspond- 
ing to weak coupling among the local displacements, the 
analytical solution is single-valued and the depinning is 
continuous. For large rj and r the analytical solution 
yields multi-valued velocity curves, reflecting the exis- 
tence of multiple sliding states, and the depinning is hys- 
teretic. As shown in the inset of Fig. 1, there is a criti- 
cal value, t]c{t), that separates single- valued from multi- 
valued solutions. The value 77 = rjc is a critical point and 
the velocity curve is expected to exhibit critical scaling. 
While the value of rjc depends on r, the existence of an 
hysteretic region at large rj, with coexistence of sliding 
and moving states and early switching (see also Fig. 1) 
occurs for all finite values of t, including t = 0. For 
T ^ 00 and 77 — !■ 00, with the ratio 11 — rj/r held fixed, 
Eq. (Q) reduces to the MFT of an overdamped elastic 
medium |^^. In this case an analytical solution is avail- 
able and the velocity vanishes linearly as F F^ 123]. 



Ti = 12, x = 4 
T| = 24, X = 8 
11 = 96, T = 32 
T] = 384,1=128 

Elastic limit, |i = 3 




FIG. 1. Analytical MFT velocity versus force curve for 
p{h) = exp(— ft) and selected values of r and rj along the 
path 77 = 3r in parameter space. As r and 77 increase along 
this path, the v{F) curve becomes hysteretic. The dashed line 
is the result for the purely elastic model with — 3, showing 
the convergence of the viscoelastic model to the elastic model 
for large r and rj. The inset shows the regions of parameters 
where depinning is hysteretic (shaded region) and where it is 
continuous (unshaded region.) 

Numerical work. We have investigated the stability of 
the branches of the analytically determined current-drive 



relationship. We performed direct numerical simulation 
of the equations of motion, for both force drive and con- 
strained mean velocity. The simulations were performed 
using two codes, for verification: a Runge-Kutta inte- 
gration and an event-driven Euler integration, with the 
"events" being crossings of a displacement Ui from one 
parabolic region to the next. The results were checked for 
insensitivity to time step At and size N. For the constant 
V constraint, the drive- velocity relationship matches the 
analytical prediction. 

In the regions where the velocity is a unique function 
of the drive, the simulation results with slowly chang- 
ing F for the force-drive curve match very closely those 
of the analytical results, which assume a constant v. In 
the presumed hysteretic region, though, the simulation 
results can be quite different. In particular, we note 
two features: mean field velocity oscillations on the lower 
branch and "early" switching, where the mean velocity 
switches from the lower to upper branch prior to the end 
of the analytically computed hysteresis region. A sam- 
ple hysteresis curve indicating early switching is shown 
in Fig. ||. We have computed the magnitude of the fluc- 
tuations in the mean velocity on the upper branch as 
a function of N: the results are numerically consistent 
with a magnitude cx N~-^/^, indicating that these fluctu- 
ations vanish as iV — > cxd. The fluctuations on the lower 
branch do not vanish in the limit of large N, however. 
These fluctuations are presumably due to an instability 
of the constant v solution in the large volume limit. We 
hypothesize, with the support of detailed analysis of the 
numerics, that nearly depinned degrees of freedom (which 
would remain pinned at constant v) are made unstable by 
velocity fluctuations and lead to an avalanche type of be- 
havior, which causes a peak in v. The magnitude of this 
instability apparently becomes large enough to drive the 
mean velocity to the upper branch before the presumed 
constant v velocity jump occurs. 

In conclusion, we have introduced a coarse-grained 
model of plastic flow that allows for slip of coherently 
pinned domains. We have solved this model analytically 
in mean field for the case of Maxwellian kernel, under the 
assumption of non-fluctuating mean velocity. We find 
that (1) the model exhibits both continuous and first or- 
der hysteretic depinning as the parameters are varied, 
(2) we can recover the case of elastic depinning in one 
limit, (3) pinned and sliding regions coexist in the hys- 
teretic regime, and (4) the mean velocity curves display 
features observed in experiments. Numerical simulations 
suggest that the behavior is much richer than suggested 
by the MF calculation and includes stick-slip-like insta- 
bilities which lead to early switching. Strong history de- 
pendence has been observed in the dc response of vortex 
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lattices in type-II superconductors [|i2|.|3| and in charge 
density waves |28 . Hysteresis in vortex lattice motion 



is most pronounced in the region of the so-called peak 
effect, where the dc response during ramp-up of the cur- 
rent proceeds via a series of jumps. These have been 
attributed to strong spatial inhomogeneities in the dis- 
tribution of vortex velocities, not unlike what is observed 
in our model We expect that in finite dimensions, 
the transition to hysteresis will be characterized by non- 
trivial universal scaling exponents | |30[ | , similar to the sit- 
uation for hysteresis in random magnets [ p2| , and that 
these exponents could be experimentally tested. 

One of us (MCM) thanks Daniel Fisher and Jennifer 
Schwarz for illuminating discussions. The work w as sup- 
por ted by N SF through grants DM R-9730678, |powre 
dmr|9805818| and |career-dmr9702242|. 
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FIG. 2. Comparison of direct numerical simulation (solid 
line) with analytic predictions (dashed line), which assume a 
constant v, for p{h) = exp(-/i), r; = 32, r = 0.8, iV = 16384, 
and a ramp rate of dF/dt — 2.5 x 10"''. The field F is cyclic 
in time. The results are in near exact agreement for much 
of the history. In the hysteretic region, on the lower branch, 
the mean field velocity occasionally spikes due to macroscopic 
events. 
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